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ABSTRACT 

Using the Kepler planet sample from Buchhave et al. and the statistical 
method clarified by Schlaufman, I show that the shorter-period super-Earths have 
a different dependence on the host star metallicity from the longer-period super- 
Earths, with the transition period being in the period range from 70 to 100 days. 

The hosts of shorter-period super-Earths are on average more metal-rich than 
those of longer-period super-Earths. The existence of such a transition period 
cannot be explained by any single theory of super-Earth formation, suggesting 
that super-Earths have formed via at least two mechanisms. 

Subject headings: methods: statistical — planetary systems — planets and satel¬ 
lites: formation — stars: statistics 


One of Kepler’s main discov eries is the existence of abundant planets wi t h sizes from 


Earth (i?®) to Neptune (4i?®) f Borncki et al.l 120 11 alibi: iBatalha et al.l 120131 iBnrke et al. 


2 OI 4 I). These so- called “super-Earths” pose challenges to the standard core accretion model 


fllda fc LinI 120041). Various models have been proposed t o explain the formation p rocess 
of such planets ( Ida fc Lin 2010 : Hansen fc Murray 2012; Chiang fc Laughlin 2013). Al¬ 
though the origin of giant planets (planetary radius i?p > 4i?®) has been reasonably well 
understood, the formation mechanism for super-Earths is still under debate. A class of 
models su g;gest that these planets are formed in situ, with, no signihcant post-assembly 
migration (iHansen fc Murray! l2 01 2l: IChiang fc LaughliiJ 120131) . Other models suggest that 
they may have a similar formation mechanism as giant planets, namely forming beyond 
the snow line where the solid material is rich and then migrating in via various interactions 
fiTerqnem fc Papaloizonll2007l: llda fc Linll2010l) . These two classes of models diverge in terms 
of the dependence on the richness of rocky material: the in situ formation mechanism gener- 
ically requires an unusually large amount of solids in the inner region of the protoplanetary 
disk, while the alternate theory does not have this restriction. Therefore, a natural way of 
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distinguishing these two theories is to look at the planet-metallicity relation of super-Earths 
as a function of orbital period. 

The planet (candidate) sample used in the present work is taken from a subgroup of 
the Kepler sar nple, which has been used to identify t he transition radius from gaseous to 
rocky planets (iRnchhave et al.ll2014l: ISchlaufmanI 120151) . The overall sample consists of 600 
exoplanet candidates orbiting a total of 405 Solar-like stars (F, G and K spectral types). 
The stellar parameters were determined based on high-resolution sp ectra gathered by the 
Kepler Follow-up Program. 0 Following the two previous studies I Rnchhave et ah 2014 : 
SchlaufmanI 120151) based on the same dataset, I also exclude small (i?p < 3i?®) planets 
that are highly irradiated (stellar flux F > 5 x 10^ J s“^ m“^) from the sample, because 
these planets could have undergone signihcant evaporatio n of their atmosph eres and thus 
no longer follow the primordial planet-metallicity relation flOwen fc Wd 120131) . This shrinks 
the sample by 23%, yielding 463 planets orbiting 324 stars . To account for t he uncertainties 
in transit depths, which are significant for small planets fIS chi an fma.nl 120151) . the error bars 
on the planetary radii of small planets {Rp < 3R^) are inflated by a factor of two with 
respect to the originally reported radii. The orbital period of each p lanet is taken from t he 
NASA Exoplanet Archive □, which co me from the Kepler pipeline f Rorucki et al.ll2011alJb : 
Ratalha et af.ll2013l: iRurke et al.ll2014j) . The sizes and orbital periods of the planets in the 
final sample are shown in Figured] as a scatter plot, color coded by the host star metallicity. 
Typical uncertainties are ~ 14% and ~ 25% on the planetary radii of large and small 
planets, respectively, and 0.08 dex on stellar metallicities; planet orbital periods are very 
well constrained (fractional error < 0.001%), so this uncertainty is not considered hereafter. 

Two populations of planets are chosen for further analysis: giant planets {Rp > 4/?®) 
and super-Earths {R^ < < 4:Rq). Planets smaller than Earth are not included in the 

smaller planet sample because studies on planet occurrence rates s uggest that these sub - 
Earth sized planets may have formed through a different channel fjPetigura et al.l 120131) : 
they are not considered as a separate population either, because only a few long-period such 
planets are detected in the current sample. 

For each planet population, the two-sample Kolmogorov-Smirnov (2KS) test is then used 
to identify whether there is a transition period. A similar method has been used to identif y 
the transition radius from gaseous to rocky planets fiRnchhave et al.ll2014j:ISchlaufmanll2015l) . 
Given an orbital period, each planet population is then divided into two groups consisting 


^See the Methods section of the original oaoer dBuchhave et al.ll2014ll for the details of the observations 
and data reduction. 


^http://exoplanetarchive.ipac.caltech.edu/. 
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of shorter-period and longer-period planets, and then the host star metallicity distributions 
of the two groups are compared, and the difference is quantihed by the 2KS test p value. 
This procedure is repeated for 10 period values from 10 to 300 days equally separated in log 
scale. To account for the uncertainties in planetary radii and stellar metallicities, I take 10^ 
realizations of the dataset, in each of which the planet radii and host star metallicities are 
randomly drawn from the Gaussian distributions dehned by the measured values and given 
error bars. The p values at each orbital period are computed for each realization, and the hnal 
reported p value for this orbital period is taken as the median, with uncertainty dehned by 
the 16 to 84 percentile values of the cumulative distribution of p values from 10^ realizations. 
As shown in the left panel of Figure [2l the p values at various orbital periods are essentially 
the same for giant planets, suggesting no distinct formation mechanisms for giant planets 
within the period range that is probed by the current sample. However, a local minimum of 
the p value distribution is found at ~ 100 days for super-Earths. Given the uncertainties in 
the p values, this local minimum does not appear at hrst sight prominent. However, these 
error bars are highly correlated, as demonstrated by Figure [3l which shows the distribution 
of differences between p values at other periods and at P = 100 d. This shows that the 
minimum near P = 100 d is indeed signihcant. Figure [3] also shows that the p value at P = 66 
day is on average equally small as that at P = 100 day, which suggests that the transition 
occurs somewhere in this range. I choose 100 day for illustration, and show the cumulative 
distributions of host star metallicities for the two groups of super-Earths separated by this 
orbital period in Figure HI For short-period (P < 100 d) and long-period (P > 100 d) 
super-Earths, the mean metallicities of host stars are 0.04 and —0.09, respectively, with 
uncertainties of 0.01 and 0.04 in these mean values. The cumulative distribution of host star 
metallicities for all giant planets is also shown in Fi gure HI with mean metall i city 0.17 and an 


uncer tainty 0.02, consistent with previous studies flFischer fc Valenti! l2005l: iRnchhave et al. 


20121 1 . 


I conduct the following tests, hrst to verify the observed local minimum in the p value 
distribution is physically real rather than the result of a random sampling effect, and then 
to investigate whether it is caused by two distinct super-Earth populations or a smooth 
transition of the metallicity-period correlation. I hrst assume a single population of super- 
Earths, with host star metallicity of each planet randomly drawn from a Gaussian distribu¬ 
tion [A/'(0.03, 0.19^), with 0.03 the mean and 0.19 the standard deviation]. Then I repeat 
the Monte Garlo process as described above to produce the p value vs. period plot, which 
is shown as blue points with error bars in the right panel of Figure [2l I then assume two 
diherent populations of super-Earths, hrst with a transition period P = 30 d and then with 
P = 100 d, the results of which are also shown in the right panel of Figure HJ In both cases 
the injected transition periods can be recovered reasonably well. The right panel of Figure [2] 
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also shows the result from the case with four different populations, with the three transition 
periods (23 d, 55 d, and 128 d) equally separated in log scale. Given the size of the planet 
sample, the four-population case represents a smooth transition on the metallicity-period 
relation, in contrast with the rapid transition for the two-population case. Compared with 
the result based on real data (the left panel of Figure [2]), the two-population model is favored 
over the smooth transition one. 


I address the influence of the systematic and selection effects before proceeding to discuss 
the theoretical implications of the results. Since Kepler planets are selected partly based 
on the signal-to-noise ratio, longer-period planets can be more easily found around brighter 
stars, which might bias on metallicities since stars with different metallicities enter the sample 
with systematically different observed fluxes. However, this selection effect turns out to be 
in the opposite direction as it would have to be in order to explain the results. In Figure [5] 
I show the scatter plot of the metallicity vs. Kepler magnitude of all super-Earth hosts, 
with colors coded by the orbital periods of surrounding planets. All targets are grouped 
into six magnitude bins, which show a tendency of gradually increasing in metallicity as 
stars become fainter. The two super-Earth populations, however, show an opposite effect: 
longer-period super-Earths tend to have hosts that are slightly fainter but more metal-poor 
than the shorter-period ones. Therefore, the inclusion of the selection effect can only make 
the two super-Earth populations more distinct. 


The results suggest that close-in super-Earths are more often found around more metal- 
rich stars, while widely separated super-Earths are more often around more metal-poor stars. 
What does it imply for the formation mechanisms of super-Earths? The in situ formation 
scenario seems to be favored over the migration scenario, since more metal-rich stars would 
presumably have more solid material at shorter distance in the protoplanetary disk and 
thus it would be easier to form super-Earths where they are now. However, the in situ 
formation scenario fails to solely explain the observed tendency for two reasons. First, if the 
in situ formation holds for all super-Earths, one would expect to see a gradually weakening 
dependence on stellar metallicity as the planet orbital period increases, rather than the very 
distinct transition period as is observed here. Another reason that also makes the in situ 
super-Earth formation very difficult is that the overall solid material required to form such 
planets is so large that it almost exceeds what is available from the protop lanetary disk could 
have at such distances f Hansen fc Murray bo 12 : Chiang fc Laughlin 2013 1. For example, the 
minimum-mass solar nebula (IHavashil 1198111 can only provide ~ d.dM® solids interior to 1 
AU. However, there are other arguments that may make the second reason unimportant, 
such as the efficiency in converting nebular material into planets being very high or most 
of these super-Earths consisting of a large amount of gas envelope. There are also theories 
that explain the super-Earths as the failed cores of giant planets. However, such theories 
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also have the problem of explaining the position of the observed transition period. 


Although a transition period is not detected in the current sample for giant planets, 
it may be that such a transition period, if it exists, is signihcantly beyond 300 days. Since 
giant planets are believed to have formed near the ice line (~ 3 AU f or Solar-like sty s) where 
plenty of solid material is accessible for building up the massive core (ILin et al.lll996h . one can 
speculate that the possible transition period may be of order thousands of days. Therefore, 
the detection of the transition period of giant planets is currently inaccessible. However, 
since the number of planets detected via radial velocity technique is growing and extending 
to longer periods, this putative transition point may become accessible in the near future. 
In addition, microlensing provides another approach to answer this question. Being sensitive 
to planets near and outside the ice line, microlensing i s probing a very dif f erent planet 


popu l ation from other t echniques flMao fc Paczvnskilll99ll: iGould fc Loeblll992l: iGould et ah 


2010l : ISnmi et al.ll201l|) . If the microlens parallax effect can be routinely measured for a 


large fraction of all planetary events, the pl anets detected via inicrolensing will have we ll 
constrained mass and distance measurements f Udalski et ahlboiS : Calchi Novati et ah 2014 1. 
When combined with the Galactic distribution of stellar metallicities, it in principle could 
also yield a planet-metallicity relation for very long-period giant planets. 
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Fig. 1.— Planet radius Rp vs. orbital period of the final planet sample; data points are 
color-coded by the host star metallicities. 
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Fig. 2.— Left panel, the p values of the two-sample Kolmogorov- Smirnov test at various 
orbital periods, for giant planets (red) and super-Earths (blue) in the iRnehhave et al.l (120141) 
planet sample. Right panel, results of the signal injection and recovery test; the blue curve is 
the result assuming a single super-Earth population. Red and black curves are the recovery 
tests of two populations, with transition period at 30 d (red) and 100 d (black); the gray 
curve shows a four-population test, which represents smooth metallicity-period relation; for 
clarity, only the blue curve shows error bars. 
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log 


p value 


p value (period = 100 d) 


Fig. 3.— The differences between p values at other orbital periods and at 100 day from 
the 10^ realizations. In large majorities of the 10® realizations, the p value at P = 100 d is 
always signihcantly smaller than p values at other periods except P = 66 d, and the p value 
at P = 66 d is on average equally small as that at P = 100 d. These results indicate that 
the local minimum is indeed signihcant, and lies near the interval P = 66 to 100 d. 
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Metallicity 


Fig. 4.— The cumulative distributions of host star metallicity for giant planets (black) and 
super-Earths with orbital period shorter than (red) and longer than (blue) 100 days. The 
average metallicities for giant planets and the two super-Earth populations are 0.17 ± 0.02, 
0.04 ± 0.01 (P < 100 d), and -0.09 ± 0.04 (P > 100 d). 
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Fig. 5.— The selection bias in the sample vs. the observed trend for super-Earths. Points 
without error bars are the positions of the super-Earth host stars in the metallicity vs. Kepler 
magnitude (brighter to fainter), color coded by the planet orbital period. This sample 
is grouped into six magnitude bins, with the vertical error bar representing the standard 
deviation of the stellar metalliciteps within each bin. The centroids for two super-Earth 
populations are marked with red (super-Earths with P < 100 d) and blue (super-Earths 
with P > 100 d). 

























